
function Cr3bp_ScnSERot(aux)
%
% 画CR3BP场景，日地旋转系
%
% 时间：2024年06月03日
% 作者：张晨
% 单位：CSU
% 邮箱：chenzhang@csu.ac.cn
%%%%%%%%%%%%%%%%%%%%%%%

% ---------------- 载入参数 ----------------
mu = aux.SB1Rot.mu;

% ----------------------- fig1 --------------------------------------
figure1;

% 太阳位置
xs = - mu;
ys = 0;
text(xs , ys , 'Sun' , 'linewidth' , 1);

% 地球位置
xe = 1 - mu;
ye = 0;
text(xe , ye , 'Earth' , 'linewidth' , 1);

% 画太阳
sunReq = 696340 / aux.dim.SB1Rot_l;
radius = sunReq;
center = [xs, ys]; % 圆心位置
theta = linspace(0, 2*pi, 100); % 构造圆的角度序列
x = center(1) + radius * cos(theta); % x坐标
y = center(2) + radius * sin(theta); % y坐标
fill(x, y, [0.9290 , 0.6940 , 0.1250]);

% 画地球
earthReq = 6378 / aux.dim.SB1Rot_l;
radius = earthReq;
center = [xe, ye]; % 圆心位置
theta = linspace(0, 2*pi, 100); % 构造圆的角度序列
x = center(1) + radius * cos(theta); % x坐标
y = center(2) + radius * sin(theta); % y坐标
fill(x, y, [0 , 0.4470 , 0.7410]);

% 画LEO
leoRmag = aux.dim.rMagLEO / aux.dim.SB1Rot_l;
plot_o([xe , ye] , leoRmag , 'k--' , 0.5);

% 画平动点
[Li_pos , ~] = Cr3bp_Li(mu);
plot(Li_pos(1,1) , Li_pos(1,2) , 'r+' , 'linewidth' , 1.2);
plot(Li_pos(2,1) , Li_pos(2,2) , 'r+' , 'linewidth' , 1.2);
plot(Li_pos(3,1) , Li_pos(3,2) , 'r+' , 'linewidth' , 1.2);
plot(Li_pos(4,1) , Li_pos(4,2) , 'r+' , 'linewidth' , 1.2);
plot(Li_pos(5,1) , Li_pos(5,2) , 'r+' , 'linewidth' , 1.2);

text(Li_pos(1,1) , Li_pos(1,2) - 0.001, 'L1' , 'linewidth' , 1.2);
text(Li_pos(2,1) , Li_pos(2,2) - 0.001, 'L2' , 'linewidth' , 1.2);
text(Li_pos(3,1) , Li_pos(3,2) - 0.001, 'L3' , 'linewidth' , 1.2);
text(Li_pos(4,1) , Li_pos(4,2) - 0.001, 'L4' , 'linewidth' , 1.2);
text(Li_pos(5,1) , Li_pos(5,2) - 0.001, 'L5' , 'linewidth' , 1.2);

xlabel('x / LU');
ylabel('y / LU');
view(0 , 90);
title('Sun-Earth Rotating')

% axisBd = 200e4 / aux.dim.SB1Rot_l;
% axis([- axisBd , axisBd , -axisBd , axisBd]);

end
